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ABSTRACT 

We present 12 mm Mopra observations of dense molecular gas towards the W28 
supernova remnant (SNR) field. The focus is on the dense molecular gas towards 
the TeV gamma-ray sources detected by the H.E.S.S. telescopes, which likely trace 
the cosmic-rays from W28 and possibly other sources in the region. Using the NH3 
inversion transitions we reveal several dense cores inside the molecular clouds, the 
majority of which coincide with high-mass star formation and Hll regions, including 
the energetic ultra-compact Hll region G5. 89-0. 39. A key exception to this is the cloud 
north east of W28, which is well-known to be disrupted as evidenced by clusters of 
1720MHz OH masers and broad CO line emission. Here we detect broad NH3, up 
to the (9,9) transition, with linewidths up to 16kms _1 . This broad NH3 emission 
spatially matches well with the TeV source HESS J1801-233 and CO emission, and 
its velocity dispersion distribution suggests external disruption from the W28 SNR 
direction. Other lines are detected, such as HC3N and HC5N, H2O masers, and many 
radio recombination lines, all of which are primarily found towards the southern high- 
mass star formation regions. These observations provide a new view onto the internal 
structures and dynamics of the dense molecular gas towards the W28 SNR field, and 
in tandem with future higher resolution TeV gamma-ray observations will offer the 
chance to probe the transport of cosmic-rays into molecular clouds. 

Key words: ISM: clouds - Hll regions - ISM: supernova remnants - cosmic rays - 
gamma-rays: observations - supernovae: individual: W28. 



1 INTRODUCTION 



(NE) and to t he sou t h (S) of the SNR (lArikawa et al.lFl 999: 



W28 (G6.4-0.1) is an old-age (> 10 4 vr: Kaspi et al.lll993h . 
mixed morphology supernova remnant (SNR) spanning 50 / x 
45 / with a distance estimated to be i n the range 1.8 to 3.3 kpc 
( e . g. I Goudis II 1 9761 : lLozinskava II 1 98 lh . The SN R exhibits non- 
thermal radio emission and t hermal X-rays (jDubner et al.l 
2000|; iRho fc Borkowskill2002h . and more recently, gamma- 

2008b|) and 
201Ch ener- 



Reach et al.l 1 2005: 



ray sources at TeV ( 10 12 eV) jAharonian et al. 
GeV (10 9 eV) (iGiuliani et aU20ld : lAbdo et al. 



gies have been discovered by H.E.S.S., AGILE, and Fermi- 
LAT telescopes respectively, pointing to high energy parti- 
cles in the region. 12 CO(1-0), 12 CO (2-1) and 12 CO (3-2) 
surveys reveal massive molecular clouds to the north east 

* E-mail: brent.nicholas@adelaide.edu.au 



^ Torres et al. ] 12003 : I Aharonian et al. 

2008bl: lFukui et al.ll2008h . Most of the CO emission appears 
centred at a local standard of rest velocity Vlsr similar to 
that inferred for W28 Vlsr ~7kms ~ *(or ~ 2 kpc) based 
on Hi studies (jVelazauez et al.ll2010h . iTorres et aU (|20Q3h 
has argued that W28 has disrupted much of this CO 
gas, giving rise to its relatively broad velocity distribu- 
tion. Notably, the NE regi on contains a rich concentration 
of 17 20 MHz OH masers (jFrail et al.lll994l : IClaussen et al.l 
1999) (with Vlsr in the r ange 5 to lSkms - 1 ), and near-IR 



rovibrational Ho emission (Reach et al. 2000; Neufeld eTaD 
120071 : Im arquez- Lugo et al. boiOh . all indicating shocked gas 



which likely results from a SNR shock interaction with 
the NE molecular cloud. The southern region contains sev- 
eral Huregions (G6. 225-0. 569, G6. 1-0.6) including the ultra- 
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2004), HESS J1714-385/CTB 37A (Aharonian et al. 2008a), 



Figu re 1. VLA 90 c m image (log scale up to 0.84 Jy/beam) from 
iBrogan et all (120061 ) of the W28 field with H.E.S.S. TeV emis- 
sion (significance contours; black dashed). Boundaries of cata- 
logued SNRs (b l ack so lid ci rcles) and 1720 MHz OH maser s from 
IClaussen et""atl (ll999T ) and iHewitt fe Yusef-Zadehl (|2009h (blue 
+'s) are also indicated. 



compact (UC)-Hn region W28 A2 (G5.89-0.39), all indi- 
cating high-mass star formation. Additional SNRs have 
also been ca talogued towards the W 28 region, namely, 
G6.67-0.42 bylYusef-Zadeh et al. I (j200d ). and G5.71-0.08 by 
Brogan et al. (2006). The ^5 arcmin resolution H.E.S.S. 
TeV gamma-ray emission is resolved into four sources. 
HESS J1800-233 is situated towards the NE region where 
a SNR shock is known to interact with a molecular cloud, 
while a group of three TeV peaks are found towards the 
south coinciding with the Hn regions (HESS J1800-240A 
and B) and the SNR candiate G5.71-0.08 (HESS J1800- 
240C). Interestingly, a recently de tected 1720 MHz OH 
maser ([Hewitt fe Yusef-Zadehl (20091 ) towards G5. 71-0.08, 
with Vlsr= 8kms"\ may also suggest HESS J1800-240C 
is tracing a SNR/molecular cloud interaction. At lower an- 
gular resolution (~5 to 20 arcmin), two Fermi-L AT GeV 
sources, 1FGL J1800.5-2359c and 1FGL J1801.3-2322c (also 
detected by the AGILE detector) appear as counterparts to 
HESS J1800-233 and HESS J1800-240B respectively. Fig- 
ure [1] compares the H.E. S.S. TeV emis s ion w ith the 90 cm 
radio (VLA) image from lBrogan et all ((2006), highlighting 
the prominent W28 SNR emission which peaks towards the 
NE interaction region, and G5. 89-0. 89 to the south. 

The TeV and GeV gamma-ray emission in the W28 re- 
gion is spatially well-matched with the molecular clouds, 
and represents the best such match outside of the cen- 
tral molecular zone (CMZ) to wards the Galactic Centre re- 
gion ((Aharonian at al. 2006). This, coupled with the old- 
age of W28, which would reduce any potential gamma-ray 
emission from accelerated electrons, suggests the gamma- 
ray emission results from collisions of cos mic-ray (CR) pro- 
tons and nuclei with th e molecular gas ((Aharonian et al.l 
l2008bl : iFuiita et al.ll2009h . W28 is a prominent member of a 
growing list of SNRs linked to TeV/GeV gamma-ray emis- 
sion spatially matched with molecular clo uds. This list in- 
cludes HESS J1745-290/SNR G359.1-0.5 ((Aharonian et al.l 



2009), 



HESS J1923 +141/SNR G49.2-0.7 dFeinstein et al 7 

and IC 443 ((Albert et al.l 120081 : lAcciari et al.ll2009h . All of 
these SNRs exhibit 1720 MHz OH masers as for W28, and 
appear to be mature SNRs (age > 10 4 yr) whereby any accel- 
erated CRs would have begun to escape into the surrounding 
interstellar medium. 

In the W28 field, the obvious source of CRs is the W28 
SNR given its prominence in many wavebands, but the other 
SNRs in the region, all with unknown distances and ages, 
may contribute to CR acceleration. Moreover, the extensive 
star formation and Hn regions associated with the southern 
molecular clouds, in particular the energetic UC-Hn region 
G5. 89-0. 39 may also contribute CRs based on recent dis- 
cussi on of protostellar particle acceleration (lAraudo et al. I 
120071 ). Some insight into where CRs are coming from can 
be obtained by looking at cloud density and emission line 
profiles in order to trace the presence and directionality of 
shocks. 



Additionally, iGabici et al.l (|2007l ) showed that the en- 
ergy and magnetic field dependent diffusion of CRs can lead 
to a hardening of the TeV gamma-ray emission spectrum 
as one looks inward towards dense molecular cloud cores 
(which typically span scale of a few arcminutes). Thus, spa- 
tial knowledge of molecular cloud density structures and 
cores towards TeV gamma-ray sources is a key step in prob- 
ing the long sought-after diffusion properties of CRs. 

Since the abundant 12 CO molecular cloud tracer, with 
critical hydrogen density ~ 10 2-3 cm -3 , rapidly becomes 
optically thick towards molecular cloud clumps and cores, 
understanding of molecular cloud density profiles and in- 
ternal dynamics can be impaired. Ideal tracers of dense gas 
such as NH3 (ammonia), CS, and HC3N are widely used due 
to their lower abundance (a factor ~ 10 _5 xCO) and higher 
critical densities ~ 10 4-5 cm -3 , which gives them a much 
lower optical thickness in dense gas. NH3 is exceptionally 
useful since a single receiver at 23-25 GHz can detect a large 
number of NH3 inversion transitions which are produced 
over a narrow bandwidth. Through satellite and hyperfine 
structure, these inversion transitions also allow the optical 
depth, and hence gas temperature and mass to be strongly 
constrained. Another desirable property of NH3 is that the 
different inversion transitions cover a wide range of excita- 
tion conditions. NH3 has therefore been detected in many 
astrophysical environments such as dense quiescent / cold gas 
and both warm to hot gas in low and high mass star forma- 
tion regions. Indeed, virtually any region containing dense 
molecular material can be studi ed with an appropriate NH3 
transition (|Ho fc Towneslll983l ). NH3 is also known to ex- 
ist in the coldest regions of molecular clouds depleting less 
rapidly from the gas phase compared to other common gas 
tracers, such as CO, whic h tends to freeze out onto dust 
grains (|Bergin et al.l l2006). 

As the next step in probing the dense cores and dy- 
namics of the W28 field molecular clouds, we have used the 
Mopra 22m single dish radio telescope in a 12 mm survey and 
single position-switched pointing covering the key inversion 
transitions of NH3 and several other 12 mm lines tracing 
high-mass star formation including H2O masers, HC3N and 
CH3OH. 
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2 MOPRA OBSERVATIONS AND DATA 
REDUCTION 

Observations were carried out on the Mopra radio telescope 
in May/ June of 2008 and April of 2009 employing the Mopra 
spectrometer (MOPS) in zoom-mode. 

Mopra is a 22 m single-dish radio telescope (31°16 / 04 // S, 
149°05 / 59 // E, 866m a.s.l.) located -450km north west of 
Sydney, Australia. The 12 mm receiver operating in the fre- 
quency range of 16-27.5 GHz, coupled with the UNSW Mo- 
pra wide-bandwidth spectrometer (MOPS), allows an in- 
stantaneous 8 GHz bandwidth. This gives Mopra the ability 
to cover most of the 12 mm band and simultaneously ob- 
serve many spectral lines. The zoom-mode of MOPS allows 
observations from up to 16 windows simultaneously, where 
each window is 137.5 MHz wide and contains 4096 channels 
in each of two polarisations. At 12 mm this gives MOPS 
an effective bandwidth of ~ 1800kms _1 with resolution 
~ 0.41 km s" 1 . Within this band the Mopr a beam FWHM 
varies from 2.4' (19 GHz) to 1.7' (27 GHz) (|Urauhart et al.l 
2010). Listed in Table [1] are some of the lines which are 
simultaneously within the bandpass in our configuration. 

On-the-fly Mapping (OTF) observations were con- 
ducted in May of 2008, and consisted of four regions, which 
are referred as A (25' x 25'), B (20' x 20') and C (15' x 15') 
to cover the TeV emiss ion peaks from HESS J1800-240 
(|Aharonian et al.ll2008bh and map D (2 0' x 20') to cover 
the Te V emission from HESS 1801-233 (|Aharonian et al.l 
2008b). We mapped each region twice, scanning once in 
right ascension and once in declination in order to reduce 
noise levels and to eliminate artificial stripes that can be 
introduced when only one scanning direction is used. It also 
allows us to check for artefacts that may occur in one scan, 
but not the other. For reference, the mapped regions are 
indicated as dashed boxes in Figure [2] 

The OTF mapping parameters we used are similar to 
those used in the the H2O Southern Galactic Plane Survey 
(HOP S), a 12 mm study of the Galactic Plane (jWalsh et al.l 
2008). Since HOPS also covered the W28 region, we have 
also included HOPS data in our mapping analysis, improv- 
ing our exposure in the mapped areas B, C & D by a factor 
of two. Map A extended beyond the Galactic latitude limit 
of HOPS (b = -0.5°) and so only partial overlap (-25%) 
exists. 

Based on mapping results and prior knowledge of 
the regions under consideration, follow-up single pointing 
position-switched deep spectra were performed in June 2008 
and in April of 2009, to provide high sensitivity to yield 
accurate measurement of the NH3 (1,1) satellite lines which 
are necessary to determine the gas temperature and density. 
As we show shortly, several regions were found to exhibit 
NH3 (1,1) and higher transitions with satellite lines appar- 
ent in the (1,1) spectra. Regions of bright NH3 (1,1) emission 
from each map were targeted in these deep spectra. In total 
there were 10 regions which were selected for follow-up spec- 
tra which consisted of 1920 s (32 min) of ON source time. 

Data were reduced using the ATNF packages 
livedata, gridzilla, ASAP and MiriadQ For mapping, 
livedata was used to perform a bandpass calibration 

1 See |http : / /www . atnf . csiro . au/ computing/ software/ for 

more information on these data reduction packages. 



Table 1. Molecular lines and the corresponding rest frequencies 
which are detectable by the MOPS spectrometer in the configura- 
tion used. The final two columns indicate whether we detect the 
line in our maps or a deep pointing. Methanol (CH3OH) masers 
are listed as type I or type II. 



Molecular Line 


Frequency 


Detected 


Detected 


Name 


(MHz) 


Map 


Deep Spectra 


H69a 


19591.11 


Yes 


Yes 


CH 3 OH(II) 


19967.396 


— 


— 


H86/3 


19978.17 


— 


Yes 


H987 


20036.32 


— 


Yes 


NH 3 (8,6) 


20719.221 


— 


— 


NH 3 (9,7) 


20735.452 


— 


— 


C 6 H 


20792.872 


— 


— 


NH 3 (7,5) 


20804.83 


— 


— 


NH 3 (11,9) 


21070.739 


— 


— 


NH 3 (4,1) 


21134.311 


— 


— 


H83/3 


22196.47 


— 


Yes 


H 2 Maser 


22235.253 


Yes 


Yes 


C 2 S 


22344.030 


— 


— 


H82/3 


23008.61 




Yes 


NH 3 (2,1) 


23098.819 




— 


CH 3 OH(II) 


23121.024 


~- 




H65a 


23404.28 


Yes 


Yes 


CH3OH 


23444.778 






NH 3 (1,1) 


23694.4709 


Yes 


Yes 


NH 3 (2,2) 


23722.6336 


Yes 


Yes 


R81/3 


23860.87 


— 


Yes 


NH 3 (3,3) 


23870.1296 


Yes 


Yes 


CH 3 OH(I) 


24928.715 


— 


— 


CH 3 OH(I) 


24933.468 


— 


— 


CH 3 OH(I) 


24934.382 


— 


— 


CH 3 OH(I) 


24959.079 


— 


— 


CH 3 OH(I) 


25018.123 






NH^ (6 6) 


25056.025 


Yes 


Yes 


CH 3 OH(I) 


25124.872 






HC 5 N(10-9) 


26626.533 


Yes 


Yes 


H89 7 


26630.71 




Yes 


H78/3 


26684.34 




Yes 


CH 3 OH(I) 


26847.205 






H62a 


26939.17 


Yes 


Yes 


HC 3 N(3-2) 


27294.078 


Yes 


Yes 


CH 3 OH(I) 


27472.501 






NH 3 (9,9) 


27477.943 




Yes 



for each row, using the preceding off scan as a refer- 
ence and applied a 1 st order polynomial fit (i.e. linear) to 
the baseline. Gridzilla re-gridded and combined all data 
from all mapping scans onto a single data cube with pix- 
els 15" x 15" x0.43 kms -1 (x,y,z). The mapping data are 
also weighted according to the relevant Tsys , Gaussian- 
smoothed (2' FWHM and 5' cut-off radius) based on the 
Mopra beam FWHM m h = 2' appropriate for the NH 3 
lines we detected, and pixel masked to remove noisy edge 
pixels. Analysis of position-switched deep pointings em- 
ployed ASAP with time- averaging, weighting by the rele- 
vant Tsys and baseline subtracted using a linear fit after 
masking of the 15 channels at each bandpass edge. In both 
mapping and position-switched data, the antenna tempera- 
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18 h 3 m 2 m l m 18 h m 59 m 58 1 
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Figure 2. NH3 (1,1) peak pixel emission [T m b ] indicating the 
four mapped regions A,B,C,D in relation to the H.E.S.S. TeV 
emission (black contours 4, 5 and 6 a levels). HOPS data are also 
included. Black/white + indicate the locations of the position- 
switched deep spectra summarised in Table [2] 

ture Tj£ (corrected for atmospheric attenuation and rearward 
loss) is converted to the main beam brightness temperature 
T m b , such that T m b = TJ r] m b where r] m b is the main beam 
efficiency. Based on the frequencies of the de tected NH 3 lines 
(^24 GHz) we assume ?7 m b = 0.6 following I Urquhart et al.l 
(|20ld ). This data reduction procedure yields an RMS er- 
ror in T m b of T rms ~ 0.05 K per channel for mapping data 
with HOPS overlap and T rms ~ 0.08 K per channel for map- 
ping data without. As a result of their increased exposure, 
posit ion- switched observations achieve a T rms of ~ 0.02 K 
per channel. 



3 RESULTS OVERVIEW 

Tabled] lists the lines detected in our mapping and position- 
switched deep spectra observations. In mapping data, the 
so called peak-pixel map for NH3 (1,1) emission is shown 
in Figure [2] along with the locations of position-switched 
deep pointings. Peak-pixel maps highlight only the bright- 
est pixel along the velocity axis (z axis) and serve as a useful 
way to search for pointlike and moderately extended fea- 
tures. 

Of the 29 molecular lines listed in Table [1] within 
the MOPS bandwidth, about half were detected. From 
the mapping observations we detected H2O, NH3 (1,1), 
NH 3 (2,2), NH 3 (3,3), NH 3 (6,6), HC 3 N(3-2), HC 5 N(10-9), 
H69a, H65a and H62a, with the criterion for detection be- 
ing a 3 T rms line peak signal. Since position-switched deep 
spectra observations were more sensitive than mapping, sev- 
eral additional lines were revealed as can be seen in Ta- 
ble [T] Maps of position- velocity (PV), integrated intensity, 
and spectra for all detected lines can be found in the online 
appendix. 

3.1 NH 3 

Presented in Figures 03 [H and [5] are integrated intensity 
and position-velocity (PV) plots. The PV plots reveal the 



Table 2. Coordinates of the Mopra 12 mm position-switched 
deep spectra observations and identification with the nearest TeV 
gamma-ray source. The position for Core 4 however, is taken from 
the mapping results since a position-switched observation was not 
taken. 



Core Name RA (J2000) Dec. (J2000) 

hh:mm:ss dd:mm:ss 

Map D - HESS J1801-233 



Core 1 




1 o.ni . cn 

18:01:59 


-16 


16 


(J4 


Core 2 




18:01:37 


-23 


26 


21 




Map A - 


HESS J1800- 


240A 






Core 3 




18:01:03 


-23 


47 


08 


Core 4 




18:02:04 


-23 


53 


04 


Core 4a 




18:01:55 


-23 


59 


05 




Map B - 


HESS J1800- 


240B 






Core 5 SW 


18:00:48 


-24 


10 


23 


Core 5 NE 




18:01:03 


-24 


08 


38 


Triple Core SE 


18:00:45 


-24 


05 


08 


Triple Core Central 


18:00:31 


-24 


04 


09 


Triple Core NW 


18:00:18 


-24 


01 


09 




Map C - 


HESS J1800- 


240C 






Core 6 




17:58:46 


-24 


09 


10 



velocity-space structure of the NH 3 (1,1), (2,2) and (3,3) 
emission regions or cores in the W28 field. A Hanning 
smoothing in velocity (width ~6kms _1 ) was applied to re- 
duce random fluctuations, and then the Galactic latitude 
axis was flattened into a single layer. In this way we show 
the intrinsic velocity location and width of the gas with- 
out confusion. Based on the PV maps, it is clear that much 
of the NH 3 emission is found in the velocity range ~-5 to 
^20kms -1 which is quite consistent with the molecular gas 
found in CO studies (|Aharonian et al.l l2008bl : iFukui et al.l 
120081 : iLiszt et al.ll2009h towards the W28 region. The four 
satellite lines of NH 3 (1,1) are clearly visible towards most of 
the cores as co- located peaks with 7 and 19kms _1 separa- 
tion from the main line for the inner and outer satellite lines 
respectively. These satellite lines spread the (1,1) emission 
over a wider -20 to 50 km s _1 range. The intensity maps also 
in Figures [3l [Hand [5] are integrated over Vlsr velocity ranges 
designed to encompass the bulk of the NH 3 emission, and 
show that it is found generally towards the TeV gamma-ray 
sources, and concentrated into clumps or cores. For simplic- 
ity we label the detected NH 3 features as Cores 1 to 6 and 
Triple Core with some of these containing sub-component 
clumps, for example, Triple Core Central, north west (NW) 
and south east (SE). Table [2] summarises the coordinates of 
the position-switched spectra towards each core, their rela- 
tion to our dedicated mapping and overlapping TeV gamma- 
ray source. 

In mapping data, most cores are seen in both NH 3 (1,1) 
and NH 3 (2,2) emission, whereas NH 3 (3,3) emission has 
been detected towards the Core 2, Core 5 SW and the 
Triple Core regions. In position-switch spectra, NH 3 (6,6) 
is detected towards Core 2 and Triple Core Central, while 
a weak detection of (9,9) is seen towards Core 2 (see online 
Figures for spectral plots). Such transitions are evidence of 
high gas temperature and potential disruption. Core 2 is 
exceptional in that it represents the well-known NE W28 
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J i i i i I i i i i I i i i i i J I 

7 6.5 6 5.5 

GLON (degrees) 

Figure 3. Top: Integrated intensity [T m b] map of Mopra NH3 (1,1) emission integrated from -20 to 50kms _1 to encompass the various 
(1,1) cores towards the TeV emission (white contours are 3, 4, 5, 6, 7 Kkms -1 ). Identified cores are labeled as described in the text. 
Dashed line indicates the Galactic Plane, dashed circles indicate the location of Cores 3 and 6 which have centre velocities outside the 
integration range. The Mopra beam (2' FWHM) is indicated in the bottom left corner. Bottom: Position- Velocity (PV) map of the peak 
pixel NH3 (1,1) emission (white contours 0.2, 0.4, and 0.6 K) indicating its velocity location and spread. For most cores the NH3 (1,1) 
satellite peaks (e.g. labeled for Core 3) may be seen as four peaks surrounding the central one along the velocity axis. 
© 0000 RAS, MNRAS 000, 000-000 
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7 6.5 6 5.5 

GLON (degrees) 

Figure 4. As for Figure [3] but showing NH3 (2,2) integrated over to 20kms _1 (white contours are 1.5, 2, 3 Kkms -1 ) and PV plot 
(white contours 0.18, 0.2, 0.3, and 0.4 K). 
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Figure 5. As for Figure [3] but showing NH3 (3,3) integrated over -5 to 20kms 1 (white contours are 1, 2, 3, 4, 5, 6, 7 Kkms 1 ) and 
PV plot (white contours 0.2, 0.4, and 0.6 K). 
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SNR shock and molecular cloud interaction region as traced 
by broad CO emission and 1720 MHz OH masers towards 
HESS J1801-233. A prominent feature here is the broadness 
and relative intensities of NH3 (3,3) and (6,6) compared to 
the (2,2) and (1,1) transitions. The Triple Core and Core 5 
comprise several resolved clumps and are linked to the en- 
ergetic Hn region G5.89-0.39 at the centre of HESS J1800- 
240B. Many of the other lines we detect in maps are found 
towards the Triple Core and Core 5, namely, bright and ex- 
tended H2O masers, radio recombination lines (RRL), and 
the cyanopolyynes HC3N and HC5N. Cores 4 and 4a ap- 
pear to trace additional Hn regions towards HESS J1800- 
240A. Core 1 is possibly a star formation site just north of 
HESS J1801-233 and the NE SNR/molecular cloud inter- 
action region. Cores 3 and 6 appear to have quite different 
velocities at ^-25kms -1 to the other cores and are likely 
not connected with the molecular gas physically associated 
with the W28 region. 

3.2 Other 12 mm Line Detections: H2O Masers, 
Radio Recombination Lines and 
Cyanopolyynes. 

A number of other 12 mm lines are detected with mapping 
data over a velocity range consistent with the W28 clouds 
and our NH3 detections. These are H2O masers, the radio 
recombination lines (RRLs) H62a, H65a and H69a, and the 
cyanopolyyne lines HC 3 N(3-2) and HC 5 N(10-9). Maps in- 
tegrated over the velocity range 5 to 20kms -1 , covering the 
bulk of the detected emission and spectra for all of these de- 
tections can be found in the online Figures. The most promi- 
nent detection of H2O masers, RRLs and cyanopolyynes is 
found towards the Triple Core region, reflecting the strong 
Hiiand star formation activity there. Since the focus of this 
paper is the dense gas traced by NH3, we only present the 
measured line parameters for the RRLs, cyanopolyynes and 
H2O masers in our online appendix. 

In the Triple Core, H2O masers are detected in both the 
SE and Central regions. The SE emission contains large ve- 
locity structure, spread over 100 km s -1 . An additional H2O 
maser is also seen towards Core 1 . According to rec ent work 
with 12.2 GHz methanol masers (jBreen et al.l l2QQ9) an evo- 
lutionary sequence for masers associated with high mass star 
formation regions has been suggested. This evolutionary se- 
quence suggests that the presence of H2O masers occurs be- 
tween ^1.5 to 4.5 x 10 4 yr after high-mass star formation, 
encompassing the onset of an initial or hypercompact (HC) 
Hn region forming ^2xl0 4 yr after the high-mass star for- 
mation. 

The radio recombination lines trace ionised gas which 
often exhibit very broad line widths as a result of likely pres- 
sure broadening and turbulence associated with ionised gas. 
The lines we see are no exception to this with several reach- 
ing line widths > 40 km s _1 extending to > 100 km s -1 in 
the Triple Core Central and SE regions. Interestingly, the 
RRL emission appears to be extended and thus may allow 
probing of the extent to which ionisation is occuring within 
the southern, central molecular cloud. As indicated in Ta- 
bleland the online appendix we detect 10 RRL transitions 
in total although only the strongest three transitions (H62a, 
H65a, H69a) are detected in the mapping data. 

The strongest cyanopolyyne HC3N and HC5N emission 



is found centred on Triple Core Central in mapping data 
but position-switched spectra reveal these molecules towards 
most of the other star formation cores (Core 1, 3, 4a, and 
6) with a possible weak detection towards Core 2. These are 
long carbon chain molecules which tend to trace the earlier 
stages of core evolution while NH3 tends to become more 
abundant at later stages ([Suzuki et al.l ll992). However re- 
cent work has indicated that HC n N(n > 3) can be produced 
under hot core conditions for short periods of time. HC3N 
is produced in large quantities at early times while HC5N is 
created and destroyed within several hundred years, m aking 
it a potential chemical clock ([Chapmann et alTll2009l ). 



3.3 CO and Infrared Comparison 

Our 12 mm line survey adds to the extensive list of molec- 
ular cloud observations devoted to the W28 SNR field. 



The Nanten telescope (Mizuno & Fukui 2004) has mapped 
W2 8 in the 12 CO (1- 0) ([Aharonian et al. I l2008bh and (2- 
1) (|Fukui et al. 1 120081) transitions, following on from ear- 
lier 12 CO (1- 0) (e.g.lArikawa et all 19991 : iDame et al.ll200ll : 



iReach et al. 1 120051: : iLiszt et al. 1 120091 ) and 13 CO (1-0) stud- 
ies ([Kim &; Kooll2003l ). More recent large-scale surveys of 
12 CO (2-1) and small scale mapping of the Triple Core re- 
gion (G5.89-3.89A and B) in 12 CO (4-3) and 12 C O (7-6) have 
been carried out by Nanten2 ([Fukui et al. .1120081 ). Our Mopra 
mapping was indeed guided by the Nanten CO results, which 
provide the most sensitive large-scale look at the molecular 
gas in the region. Figure[6]compares the Nanten2 12 CO (2-1) 
image with the NH3 (1,1) and (3,3) emission from our Mopra 
observations, No beyama 12 CO ( 1 -0) an d JCMT 12 CO (3-2) 
observations by Arika wa et al. (1999) as well as the TeV 
gamma-ray emission with H.E.S.S. The 12 CO (2-1) emission 
spatially matches well the brightest three TeV gamma-ray 
peak s as highlighted previou sly for the 12 CO (1-0) emission 
by lAharonian et al.l (|2008bl ) . This match is quite striking 
towards the Core 2/HESS J1801-233 region where broad 

12 CO (3-2) or shocked/disrupted gas tends to lie inward in 
the direction of W28 compared to the more quiescent and 
relatively narrower 12 CO(1-0). Here, the NH3 (3,3) emis- 
sion reveals for the first time the dense and disrupted core 
of the shock-compressed NE molecular cloud. Towards the 
Triple Core region the Nanten2 12 CO (2-1) emission is re- 
solved into two peaks assocated with the Hn regions G5.89- 
3. 89 A and B respectively. Core 5 NW and SW are also traced 
by 12 CO (2- 1) enhancements, w hich are also clearly seen in 

13 CO (1-0) ([Kim fc Kooll2003l ). Our Mopra NH 3 emission 
also appears to resolve peaks associated with G5.89-3.89A 
and B, and Cores 5 NE and SW. 

Figure [7] presents the infrared (IR) Spitzer GLIMPSE 
and MIPSGAL three-colour (RGB=24/8/3.6 /mi) image of 
the W28 region with Mopra NH 3 (1,1) and H.E.S.S. TeV 
gamma-ray contours. These IR bands are tracers of poly- 
cyclic aromatic hydrocarbons (PAHs) and dust emission, re- 
vealing the complexity of the W28 region in hosting several 
star formation and Hn regions in quite likely several different 
evolutionary stages. Table [3] summarises the likely counter- 
parts identifed with our detected cores. 

Further discussion of the individual cores and compar- 
ison of previous molecular line studies toward them can be 
found later in §5\ 
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Figure 6. Left: N anten 12 CO (2-1) image [Kkms -1 ] 
(|Fukui et al. 1 120081) in log scale, with contours of Mo- 
pra NH 3 (1,1) (white) and H.E.S.S. TeV gamma-ray 
significance (black-das h ed). 1720 MHz OH mase r s from 
IClaussen et aD (|l999l ): iHewitt fe Yusef-Zadehl (j2009T ) 
are also indicated (blue/white +). Above: Core 2 zoom 
in linear scale with contours of 12 CO ( 1-0) (magenta- 
dashe d) and 12 CO (3-2) (blue) from [Arikawa et all 
(Il999h. SNR diameters for W28 and othe r SNRs 
(|Yusef-Zadeh et al.lliuOol : iBrogan et al.ll2006h are in- 
dicated in both panels. 
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Table 3. Likely counterparts to the various cores detected. Coun- 
terparts have been searched for within a 2' radius using the 
SIMBAD^ astronomical database. Where a star is listed the spec- 
tral type is also given. 



Core Name 



Counterpart Name 



Figure 7. Spitzer GLIMPSE/MIPSGAL three colour 
(RGB=24/8/3.6 ^m; MJy/sr in log scale) image with con- 
tours from Mopra NH3 (1,1) emission (linear white contours 
> 1.5K) and H.E.S.S. TeV gamma-ray significance (yellow 
dashed). 



Core 1 
Core 2 
Core 3 
Core 4 
Core 4a 
Core 5 SW 
Core 5 NE 

Triple Core SE 
Triple Core Central 
Triple Core NW 
Core 6 



IRAS 17589-2312 1 

W28 SNR/Mol. cloud interaction 

HD 313632 (B8IV) 

Hn 6.225-0.569 2 

Hn G6.1-0.6 3 / IRAS 17588-2358 

IRAS 17578-2409 

HD 164194 (B3II/III C) 

IRAS 17578-2409° 

Hn G5.89-0.39A 4 / W28-A2 

Hn G5.89-0.39B 4 

V5561 Sgr (M) / IRC-20411 5 

IRAS 17555-2408 6 



Bronfman et al. (1996) 

Lockman (1989) 

Kuchar & Clark ( 1997) 
4. see Kim 8z Koo (2001) and references therein. 
5. 1 Johnst on et al. ( 19731) . 
| http : //simbad . u- strasbg . f r/simbad/ 

a. 2.8 7 distant. 

b. 2.3 ; distant. 
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4 ANALYSIS OF NH 3 EMISSION 

We describe here procedures used to estimate the gas pa- 
rameters such as optical depth, temperature, mass, density, 
and related dynamical information towards the NH3 cores. 
Table^l] (with statistical errors in the online appendix) and 
[5] summarise these results utilising the NH3 (1,1) and (2,2) 
spectra for an analysis assuming various core sizes, and ad- 
ditional results treating Core 2, Triple Core and Core 5 as 
extended clouds well beyond the size of the 2' FWHM beam. 
We also apply a detailed radiative transfer model to Core 2, 
given its apparent high gas temperature and non-thermal 
energy. 



4.1 Line Widths 

The FWHM of the NH3 main line, Aui/2, is a useful mea- 
sure of the total energy associated with the core or clump. 
Broader lines will result from regions with higher temper- 
atures or some additional dynamics. The purely Maxwell- 
Boltzmann thermal line width FWHM A vth expected from 
a gas at temperature T is given by 



A^th 



ln(2) k T _ x 

— [km s J 



77iNH 3 



(1) 



where k is Boltzmann's constant and mNH 3 is the mass 
of the NH3 molecule. For example, a thermal line width 
of ^0.16 kins" 1 is obtained for a temperature of 10 K, 
as might be expect ed in typically cold dense NH3 cores 
(|Ho fe Towneslll983l ). The line FWHM, Av 1/2 , of each core 
was estimated from a Gaussian fit to the central or main 
peak of the emission with additional Gaussians to fit each 
of the four satellite lines, which are generally resolved in the 
(1,1) transition. This five- Gaussian fit function can be seen 
applied to the Core 1 and Core 2 spectra in Figure [8] The re- 
sults in Table[4]show that A v 1 / 2 for all cores is considerably 
wider than that expected from purely thermal broadening, 
suggesting additional non-thermal or kinetic energy which 
dominates over broadening from the instrumental response, 
which is considered negligible. 

4.2 Gas Parameters 

The (1,1) satellite lines were clearly resolved in all cores 
except in Core 2 which is instrinsically very broad, lead- 
ing to blending of the main and satellite lines (see Fig- 
ure [8]) . Using the relative brightness temperatures of the 
main and satellite peaks, the optical depth can be derived 
for each (J, K) in version transit i on by numerically solv- 
ing Equation 2 of iBarrett et al.l (| 19771 ). A weighted aver- 
age f(J,K,m) (with weight given as 1 /peak-error accord- 
ing to the Gaussian fit) of the four satellite-derived optical 
depths is then calculated, and the total optical depth t]°k 
of the transition is estimated by accounting for the frac- 
tion of intensity in the main line (where /( 1,1) =0.502, and 
/(2,2)=0.796). Without needing to resolve the (2,2) satel- 
lite lines, the main line optical depth of the (1,1) transi- 
tion can be used to infer the main line optical depth of 
the (2 ,2) transition using the method of lUngerechts et al.l 
(1986]). As for the (1,1) transition, a weighted average may 
be used to derive the total (2,2) optical depth t%°2- The 
line FWHM is of critical consideration when determining 



NH 3 (6,6) 
NH 3 (3,3) 



2.5 



NH 3 (9,9) 
NH 3 (6,6) 




Figure 8. NH3 (n,n) position-switched spectra for Core 1 (Left) 
and Core 2 (Right). Gaussian fits to the (1,1) and (2,2) spectra 
used to estimate gas parameters (see Table [4|). The (1,1) fits are 
shown as red solid lines and the (n,n) spectra for n^2 are offset 
by a constant value for clarity. 



the temperature of the gas, as a key assumption in this step 
is that both the (1,1) and (2,2) emission probe the same 
volume of gas. The similarity of FWHM for the (1,1) and 
(2,2) transitions in Table [4] does however suggest that both 
transitions originate from the same general volume of gas, 
allowing the r otational tempera t ure to be calculated using 
the method of lUngerechts et al~l (|l986h . Generally the rota- 
tional te mperature is an underestimate of the kinetic tern- 

Ijasi), 

/ onn/i 



perature ( Ho & Townes 



pression of lTafalla et al. 



however, the analytical ex- 
(|2004 pg. 211) has been used to es- 
timate Tk from T r ot • This is believed to be accurate to within 
5% of the real value of Tk for temperatures in the 5-20 K 
range. Importantly, this method for obtaining temperatures 
is only considered valid for Tk up to ^40 K. To obtain the 
column density for the NH3 e mission, we use the (1, 1 ) tran - 



sition and follow Equation 9 of I Goldsmith &; Langerl (|1999 
Since the (1,1) represents only a fraction of the NH3 gas, the 
total NH3 column density Anh 3 can be obtained from the 
(1,1) column density and the partition function Q(T). In 
this work we assume contributions up to the 5 th term in the 
series expansion. 

To determine the total mass M and molecular hydro- 
gen number density 71h 2 (hereafter referred as just den- 
sity) of the various cores we need to assume the core in- 
trinsic radius R and the abundance ratio Xnh 3 of NH3 to 
H2. We note that the abundance ratio is known to depend 
on the gas temperature T. Rati os quoted in literature ex- 
tend from ~1(T 7 , few xlO" 8 (IWalmslev fe Schilkel \l9s3i : 
lOtt et al.ll2005h . down to - 10 10 lOtt et al.l (|2010h in the 
Large Magellanic Cloud. For typical infrared dark clouds 
(IRDCs) with T < 20 K an average value of - 4 x 10~ 8 
is found ([Pillai et al.l l2006). However with the exception of 
Core 2 and Core 6, all NH3 cores in the W28 field are identi- 
fied with star formation IR sources as indicated in Figure [71 
and Table [3] In hot cores, with (T > 100 K), abundance 
ratios of 10~ 6 or even as high a s 10 -5 have been suggested 
(|Walmslev fe Schilkell983l . ll992l and references therein). We 
adopt the value 



Xnh 3 = AW 3 / N H2 = 2 X 10" 



(2) 
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following Istahler fe PalTal (|20Q5 . Table 5.1) as the temper- 
atures of our detected cores are only slightly higher than 
those typical of IRDCs. The total hydrogen mass of the core 
is then estimated by factoring in the projected area of the 
emitting region: 

fK 2m H 7ri? c 2 m AW 3 /(XNH3 M ) [M ] 



M 



(3) 



for the core intrinsic radius expressed in cm units R cm 
and 7 71h the hydrogen mass (kg). Following lUrquhart et al.l 
(|2010h we apply the beam dilution factor /, given by 1/ f — 
In 2 (2i?/# m b) 2 , and beam coupling coefficient K — x 2 /[l — 
exp(— x 2 )] for x — y/T/f in order to account for various 
source radii. Therefore the product fK = 1/[1 — exp( — 1//)] 
with fK — )• 1.0 for increasing R. 

The H 2 molecular number density nn 2 follows by ac- 
counting for the emitting volume of emission, assumed to 
be a sphere with radius R C m in cm. 



nH 2 MM /(2ra H 4/37ri?cm) [cm 



(4) 



Except for Core 2 and possibly the Triple Core and Core 5 
complexes, the NH3 emission appears to be point-like and 
thus has an intrinsic diameter upper limit approaching the 
2' Mopra beam FWHM, or R mh = 0.56> mb = 0.6 pc radius 
assuming a 2 kpc distance (consistent with the W28 SNR). 
Although R = 0.1 pc is t ypically noted as th e intrinsic ra- 
dius for cold dense cores Ho fe Townes II 19831 ). we assume a 
source radius of 0.2 pc as a default. The estimates for mass 
and density assuming a point-like source radius of R = 0.2 pc 
are listed in Table [5] In order to allow for varying core radii, 
we also include scaling factors for the mass and density re- 
sulting from the source radius R dependence in fK, source 
projected area, and volume. 

In considering Core 2, Triple Core and Core 5 as ex- 
tended regions, we calculate gas parameters from their 
NH3 (1,1) and (2,2) spectra averaged over elliptical regions 
and for pixels T m b ^ 0.18 K. The extended mass is esti- 
mated as in Eq. [3] using the column density averaged over 
the extended region of radius R C m in cm, but with an extra 
term ?7mb/?7xb = 0.86 to account for the extended Mopra 
beam efficiency rj^h = 0.7 appropriate for the NH3 lines 
([Urquhart et al.lfeoid ). Density then follows from Eq. [Has 
before. Results, including dimensions of the elliptical regions 
are given in Tables [H and [5] Additionally for Core 2, we use 
the NH3 (1,1) to (6,6) spectra (position-switched and map- 
ping) in more detailed radiative transfer modeling discussed 
in M4.4I to estimate gas parameters. 

Under the assumption that the cores are in gravitational 
equilibrium with their thermal energy, their pure molecular 
hydrogen virial masses M wlv may also be estimated: 



M vir = kR(Av 1/2 



[M G 



(5) 



for the source radius R (pc) as before and Avi/ 2 the line 
FWHM (kms -1 ). The factor k depends on the assumed 
mass density profile of t he core p(r) wit h ra dius r. For a 
Gaussian density profile IProtheroe et al.l (|2008h calculates 
k = 444, in contrast to other situations s uch as a constan t 
density (k = 210) and p oc r~ 2 (k = 126) (|MacLaren II 19881 ). 
Although a Gaussian profile is quite likely, we quote here in 
Table [5] the virial masses bounded by the Gaussian and r~ 2 
density profiles with a source radius i?=0.2pc. Additionally, 
an overestimate of the true core line FWHM can result for 
optically thick lines. Given the similarity of the NH3 (1,1) 



and (2,2) linewidths, and that the (2,2) optical depth is gen- 
erally less than unity (see Table SJ), we use the (2,2) FWHM 
in the virial mass calculation. The exception is for Core 2 
analysis in which case we use the (3,3) linewidth. Overall, 
the virial masses we derive here could be considered upper 
limits, especially in the case of a Gaussian density profile. 
Nevertheless, the virial mass serves as an important guide 
in understanding the stability of the cores. 

By far the dominant systematic error in our mass and 
density estimates arises from the uncertainty in the NH3 
abundance ratio Xnh 3 - Given the range of ratios quoted in 
literature for cores of similar temperature to ours, we quote 
systematic errors of a factor 2 to 5 for Xnh 3 , which feed 
directly into mass and density. Given that most of the core 
masses we derive are in agreement with their virial mass 
range (Table O, our choice of Xnh 3 = 2 x 10 -8 appears 
reasonable. 



4.3 Velocity Dispersion: Core 2 and 
Triple Core/Core 5 

The detection of broad NH3 emission primarily towards 
Core 2 and Triple Core suggests active disruption of the 
molecular material. The dynamics of the NH3 gas can be 
probed by looking at the velocity dispersion across a cloud 
core at each pixel v rm s weighted by its NH3 intensity (T m b ) , 
in addition to the position-velocity information (see online 
appendix). For each pixel with intensity above a reasonable 
threshold, in this case 0.18 K Or ^3.5 Trms 5 

the velocity dis- 
persion is calculated as: 



_ // T mh (v) (v - v) 2 dv 
~ V J T mh (v)dv 



[km e 



(6) 



for v = j vT m h (v)dv/ j T m b (v)dv the intensity- weighted 
velocity. Results are presented in Figure [9] for the Core 2 and 
the Triple Core/Core 5 regions using the NH 3 (1,1), (2,2) 
and (3,3) transitions. In the (2,2) and (3,3) transitions a 
wide -50 to 50 km s _1 Vlsr velocity range encompassing the 
bulk of the emission was considered. For the (1,1) emission 
however, the strong satellite lines can contaminate Eq. [6] 
and thus a restricted Vlsr range was used for Core 2 (5 
to 15 kms -1 ) and Triple Core/Core 5 (5 to 20 kms -1 ). To 
further remove fluctuations, data cubes have been Hanning 
smoothed with velocity width ^2 kms" 1 . For the Core 2 re - 
gion, the 1720MHz OH masers from IClaussen et al.l (|l999h 
are indicated in Figure [9] in order to outline the regions 
where the W28 SNR shock is interacting directly with the 
NE molecular cloud. For the Triple Core/Core 5 region, the 
Hn regions G5.89-0.39A and B are indicated. 

While the magnitude of the dispersion varies between 
the three NH3 lines in both cores, the dispersion maps indi- 
cate the contrast between Core 2 and the Triple Core/Core 5 
region. Core 2 has the most velocity dispersion in the (3,3) 
line, indicated by both the peak magnitude and the physi- 
cal area, both of which get progressively smaller in the (2,2) 
and (1,1) lines. However, the Triple Core/ Core 5 region has 
the most dispersion in the (1,1) line, and progressively less 
in the (2,2) and (3,3) lines. Interestingly, the spatial loca- 
tion of the peak of the disruption in Core 2 moves radially 
outward in the direction of the W28 SNR shock, from the 
western side of the core (in the (3,3) line) toward the eastern 
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Table 4. NH3 gas parameters derived from NH3 (1,1) and (2,2) spectra. Results for point-like analysis are taken from a single position 
(either from position-switched data or from mapping data in the case of Core 4). Extended source parameters are taken from spectra 
averaged over an elliptical region given below and for pixels satisfying a T m b ^ 0.18 K masking level in mapping observations. Columns 
from left to right are core name, integrated intensity for (1,1), (2,2) and (3,3), rotational temperature T ro t , kinetic temperature Ty. , 
NH3 (1,1, m) main line velocity Vlsr 5 main line FWHM Av 1 / 2 for (1,1), (2,2), and (3,3), total column density A^ NH3 , and total optical 
depth rl°i and t\°2 f° r (1,1) and (2,2) respectively. Statistical errors are given in the online appendix. 
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3.8/1.0/0.8 


20.1 


25.0 


+9 


4.5/4.7/5.2 


15.0 


1.2/0.7 


f 5% systematic errors also apply. 














% I T mb dv and Avj/ 2 for 


NH 3 (l,l)/(2,2)/(3,3)/(6, 


6)/(9,9) 












1. For ellipse 5.2pcX3.5pc diam. (spherical radius 2.1pc); pos. 


angle +30° 


; RA 18:01:41 Dec 


-23:25:06 






2. For ellipse 7.0pcX2.4pc d 


iam. (spherical radius 2.1pc); pos. 


angle +30° 


; RA 18:00:49 Dec 


-24:10:23 






3. For ellipse 6.6pcx3.1pc d 


iam. (spherical radius 3.2pc); pos. 


angle -30° ; 


RA 18:00:30 Dec - 


24:03:09 







side of the core (in the (2,2) and (1,1) lines). For the Triple 
Core/Core 5 region, the velocity dispersion is always peaked 
toward the centre of the cores, with no evidence for disrup- 
tion form an external source. This evidence would suggest 
that the W28 SNR shock has disrupted Core 2 but not yet 
reached the Triple Core/Core 5 region. 

4.4 Radiative Transfer Modeling of Core 2 

The broad NH3 (3,3) line with Av 1 / 2 ^13kms -1 and its 
high line strength relative to (1,1) and (2,2) from Core 2 
are somewhat beyond those expected from a purely thermal 
distribution as indicated in £14.11 Such broad line profiles 
suggest a large amount of energy, in this case non-thermal 
energy from the SNR shock, has been deposited into the 
cloud. The velocity dispersion NH3 (3,3) image for Core 2 
(Figure [9]) and the many 1720 MHz OH masers in the region 
indicating shocked gas support the notion that the broad 
NH3 is also the result of SNR shock disruption. The analysis 
outlined earlier in § 14.21 for the (1,1) and (2,2) emission as- 
sumes that it comes from quiescent, cold or cool, dense cores. 
This assumption clearly doesn't hold for the Core 2 region, 
and as such the mass and density estimates for Core 2 in 
Table [5] are likely underestimates which would pertain pri- 
marily to the cooler gas component. We therefore turn to 
more detailed radiative transfer modeling of the position- 
switched (1,1), (2,2), (3,3), and (6,6) spectra taken towards 
the peak of NH3 (3,3) emission (Figure [8]) and spectra av- 



eraged over the Core 2 region for a first look at the gas 
parameters towards Core 2. 

The radiative transfer code, MOLLIE0, can deal with 
arbitrary 3D geometries, but as a first step in obtaining 
the indicative properties of Core 2 in this paper, we mod- 
eled the emission as arising from a sphere with a constant 
temperature, density and non-thermal velocity component, 
taking the NH3 spectra corrected for the Mopra aperture 
main beam efficiency of r] m b — 0.6. The NH3 to H2 abun- 
dance ratio was fixed at 2 x 10 -8 as in our earlier anal- 
yses, and a source radius of 2.1 pc (distance 2kpc) was 
chosen based on the extent of the NH3 (3,3) intensity af- 
ter a T m b ^ 0.18 K cut. Models were constructed with H2 
densities, temperatures and non-thermal linewidths ranging 
from 10 3 to 10 8 cm" 3 , 10 to 400 K and 0.5 to 40kms"\ 
respectively. Radiative transfer modeling was then used to 
generate synthetic data cubes with a velocity resolution of 
O.lkms- 1 for the NH 3 (1,1) , (2,2), (3,3) and (6,6) emis- 
sion. These were then convolved with 2D Gaussian pro- 
files at a spatial scale corresponding to the 2' FWHM of 
Mopra. The synthetic spectra at each transition were fit 
to the observed spectra (weighted by the signal-to-noise 
of each transition) and reduced-x 2 values returned for the 

2 See | Ketol fl990) for a description of the MOLLIE code, 
Keto et al.l (l2004f) f or a description of the li ne-fitting and 
Keto fc Zhangl (|2010h and iLongmore et al.l (|201(t ) for recent ex- 
amples of work using the code 
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Figure 9. Intensity- weighted velocity dispersion u rms [kms -1 ] of the NH 3 (3,3), (2,2) and (1,1) emission. The (3,3) and (2,2) emission 
is considered over the -50 to 50 kms -1 VLSR ran g e > whilst the (1,1) is considered over the 5 to 15 kms -1 (Core 2) and 5 to 20 kms -1 
(Triple Core/Core 5) ranges to avoid the strong satellite lines contaminating the i; rms calculation. Pixel channels are masked below a 
T m b value of 0.18 K (~3.5T rm s) threshold. Dashed black contours indi cate the H.E.S.S. TeV gamma-ray significance. For Core 2 the 
blue/white +'s indica te the positions of t he 1720 MHz OH masers from IClaussen et al.l £l999), and red x's indicate positions of broad 
line regions denned bv lReach et al.l (|2Q05h . The large and small solid circles represent the radio boundaries of the SNRs W28 and G6.67- 
0.42 respectively. For Triple Core/Core 5, blue/white +'s indicate G5.89-0.39A (Hn ) and B (UC-Hn). The Mopra beam 2' FWHM 
is indicated on the Core 2 NH3 (2,2) image and applies to all other images shown here. The magenta solid ellipses define regions for 
extended source mass and density estimates in Table The implications of the more central dispersion seen in the Triple Core versus 
the generally asymmetric dispersion seen in Core 2 are discussed in the text. 



goodness-of-fit. Simulated annealing with 10,000 models was 
used to search through the 3D parameter space to min- 
imise x 2 an d find the best-fit model. This method is inher- 
ently robust against becoming trapped in local, rather than 
global minima in parameter space. However, to determine 
the robustness of the best-fit model we ran the fitting 20 
times with widely separated initial start values and incre- 
ments. For the position-switched spectra, the best-fit model 
yielded a hydrogen atom number density, temperature and 
non-thermal linewdith of 10 3,45 cm -3 , 95 K and 7.4kms _1 , 
respectively, giving a Core 2 mass of —2700 M©. For the 
mapping- averaged spectra, the best-fit model yielded a den- 
sity, temperature and non-thermal linewdith of 10 312 cm -3 , 
60 K and 7.5 kms -1 , respectively, giving a Core 2 mass of 
-1300 M . 

Figure [10] shows the NH 3 (1,1), (2,2), (3,3) and 
(6,6) Mopra spectra from position-switched and averaged- 



mapping data overlayed with the synthetic spectra from the 
best-fit model. The weakly detected (9,9) spectrum was not 
included in both cases. Considering the simplicity of the 
model, the synthetic spectra match the position-switched 
data well. Differences between the model and data provide 
insight into the underlying source structure of Core 2. The 
single non-thermal velocity contribution to the linewidth 
works well for the higher transitions, but can not account 
for the narrow linewidth component of the NH3 (1,1) emis- 
sion from both sets of spectra used. Similarly, assuming a 
single temperature for Core 2 underestimates the NH3 (6,6) 
emission, especially for the mapping- averaged spectra, which 
yield a lower density and temperature compared to position- 
switched results as expected. Given these issues, the Core 2 
extended masses derived here are likely to be underesti- 
mates and we conservatively quote the lower of the two 
densities and masses in Table [5] More detailed modeling 



© 0000 RAS, MNRAS 000, 000-000 



14 B. Nicholas et al 



NH3(1,1) Position-Switched 




NH3(2,2) 




NH3(3,3) 




o NH3(6,6) 




-50 50 

Velocity [km/s] 
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Figure 10. MOLLIE radiative transfer modeling fits to the NH3 
spectra from Core 2 outlined in £14.41 Left-Column: Fits to 
position-switched deep spectra towards the peak NH3 (3,3) emis- 
sion. Right-Column: Fits to mapping spectra averaged over pix- 
els with T m b ^ 0.18 K. In all panels the best-fit model is overlaid 
as the red solid line. 



to simultaneously fit the cooler gas traced by the narrow 
linewidth NH3 (1,1) emission and the hotter gas traced by 
the NH3 (6,6) emission will be the focus of later work. Errors 
in the fitted quantities are not yet estimated due to the of- 
ten imperfect fits to each spectra, and will also be discussed 
in later work. 



5 DETAILED DISCUSSION OF CORES 
Core 1: 

Situated at the northern boundary of HESS J1801- 
233 in the vicinity of the giant Hil region M20. Core 1 
exhibits NH 3 (1,1), (2,2) and (3,3) emission centred at 
VLSR=+21kms -1 with T k ~ 18 K. The (1,1) and (2,2) 
lines are detected in mapping while the (3,3) line is only 
seen in the deep pointing. CS (2 -1) was earlier detecte d 
at a similar Vlsr value found by iBronfman et al. I Jl996), 
who first suggested the link to the IR source IRAS 17589- 
2312. Based on the FIR colour ratio, these authors have 
also sugge sted IRAS 175 8 9-231 2 is indicative of an UC- 
Hn region. iFaiindez et al. I (|2004r ) estimate a core mass of 
M = 210 Mq, radius ~0.25pc (for 3.8 kpc distance), and 
density nn 2 = 3.4 x 10 4 cm -3 based on their SIMBA 1.2 mm 
continuum mapping. Assuming a 0.25 pc radius, the mass 
and density estimates from our observations (500 Mq and 
1 5.0x 10 4 cm~ 3 ) are a f actor 2 to 4 times higher than those 
of lFaundez et al.l (|2004h , but in general agreement given the 
systematic uncertainties arising from the NH3 abundance 



ratio. A similar core radius is also quoted by iLefioch et al.l 
(2008) who studied this core using a variety of CO, CS, and 
HCO + transitions. They suggest that star formation activity 
in this core is very recent which could have been triggered 
by the SNR W28. The H2O maser seen tow ards this core 
has also been discussed ([Codella et al.lll995h . 

Core 2: 

As already described, Core 2 spatially coincides very 
well with the TeV gamma-ray sou rce HESS J1801-233 
and broad- line CO l ine observations (lArikawa et al.l 1999; 



Torres et al. 2003 



IReach et al.l |2005k lAharonian et al.l 



2008bl: lFukui et alJIiood ). This molecular cloud is known to 



be shock-disrupted as e videnced by the p r esence of many 
1720 MHz OH masers ([Frail et al.l Il994l : IClaussen et~ah1 
1999). Further detailed studies of the shocked gas have been 
carri ed out in CO, CS and IR H2 lines by IReach et al.l 
(2005). Our Mopra observations reveal NH3 inversion tran- 
sitions from (1,1) up to (9,9) with very broad linewidths. 
The strongest emission with j T m b dv = 7.4Kkms _1 is 
found in the broad A v 1 / 2 = 12.8 kins" 1 (3,3) transition, 
which bears a close resemblance to the CO peaks. The ve- 
locity dispersion image of the (3,3) line in Figure [9] shows 
clearly that the cloud disruption originates fr om the western 
or W2 8 SNR side, supporting the results of lArikawa et al.l 
(1999) who showed that the broad 12 CO (3-2) emission is 
found preferentially towards the W28 side in contrast to 
12 CO (1-0) which extends radially further away. The Nan- 
ten2 12 CO (2-1) emission likely traces a mixture of shocked 
and unshocked gas and detailed velocity dispersion studies 
of this (2-1) emission are currently underway. The lack of 
strong broad-band IR features towards Core 2 (Figure [7]) 
suggests that the NH3 excitation is not due to star forma- 
tion processes, although some fraction of the weak IR emis- 
sion seen here is due to s hocked H2 (jNeufeld et al. 1 120071 ; 
iMarquez-Lugo et al. , 2010). We also note there is a weak 
HC3N feature seen in the deep pointing which is an indica- 
tor of hot gas-phase chemistry. It is also quite striking that 
a grouping of the OH masers appear to surround the re- 
gion where the (3,3) emission is most disrupted (Figure [9}, 
radially away from the W28 SNR. This would be further evi- 
dence in support of the W28 SNR as the source of molecular 
cloud disruption, and overall would tend to disfavour physi- 
cal influence from the neighbouring SNR G6. 67-0. 42, which 
at present has an unknow n distance. The b road molecular 
line regions discussed by IReach et al. (2005) also generally 
cluster towards the broadest NH3 emission (see Figure [9}. 

The broadening of the NH3 (3,3), (6,6) and probably 
(9,9) lines are dominated by non-thermal component (s) and 
we can estimate the additional kinetic energy Wkin required 
to achieve this from: 



W kin = l/2M(A^ kin ) 2 



where M is the mass of the broad-line gas and At^m is the 
FWHM [ms _1 ] of the line due to additional non-thermal 
kinetic processes. Using the non-thermal line width A^kin = 
7.5kms _1 and mass lower limit M = 1300 Mq from our ra- 
diative transfer modeling in £14.41 of mapping- averaged NH3 
spectra we therefore calculate Wke > 0.7 x 10 48 erg. This 
energy lower limit is within a factor few of the kinetic energy 
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Table 5. Mass M and molecular hydrogen number density nn 2 estimates (Eqs. [3]and[4j) for the various cores from JVnh 3 derived using 
NH3 (1,1) and (2,2) parameters of Table 0] Point source analysis assumes a spherical emission volume of radius R = 0.2 pc (with the 
effect of source radius scaling indicated) whilst for extended source analysis spectra are averaged over an elliptical region (see Table [4j) 
using pixels satisfying T mb ^ 0.18 K. The included range of virial masses M v [ r is bounded by radial density profiles following r~ 2 and 
Gaussian laws and assumes the NH3 (2,2) linewidth Av 1 / 2 (except for Core 2 which assumes the (3,3) linewidth). The radiative transfer 
results for Core 2 using the MOLLIE code applied to the NH3 (1,1) to (6,6) spectra are also given. 



Core Name 


M 


7lH 2 


M vir 




[M ] 


[10 3 cm- 3 ] 


[M ] 


Point Source Analysis 


(R = 0.2 pc) 




Core 1 


495 


287.9 


220 - 760 


Core 2 


240 


140.8 


4200 - 14700 f 


Core 3 


640 


370.2 


280 - 980 


Core 4 


190 


112.3 


80 - 290 


Core 4a 


75 


43.5 


230 - 810 


Core 5 SW 


350 


202.5 


140 - 480 


Core 5 NE 


335 


196.1 


100 - 360 


Triple Core SE 


475 


276.8 


330 - 1170 


Triple Core Central 


320 


185.1 


520 - 1820 


Triple Core NW 


130 


76.7 


110 - 400 


Core 6 


910 


530.0 


220 - 760 


Point source mass/density scaling vs. radius R (pc) 




R (pc) 


0.10 0.15 0.20 0.25 0.30 0.35 0.40 


Mass M(R) /M(0.2pc) 


0.97 0.98 1.00 1.02 1.05 1.08 1.12 


Density n H2 (R) /n H2 (0.2pc) 7.77 2.33 1.00 0.52 0.31 0.20 0.14 




Extended Source Analysis 




Core 2 


1600 


0.8 


44400 - 152600 f 


Core 2 (MOLLIE) 


>1300 


>0.7 


14500 - 51200 x 


Core 5 


1300 


0.7 


4000 - 14200 


Triple Core 


3300 


0.5 


8900 - 31400 



f Using from the (3,3) emissk 



l Table [4] 



y l/2 



as the non-thermal line width from MOLLIE. 



(~ 3 X 10 48 erg) deposited i nto the 2000 Mq of gas traced 
by shocked 12 CO (3-2) from lArikawa et al.lll999l . 

Over the to 12 km s _1 Vlsr range, for which the Nan- 
ten 12 CO (1-0) emission shows excellent overlap with the 
TeV gamma-ray source HESS J1801-233, the mass of the 
NE molecular cloud is ~ 2 x 1O 4 M . Over a wider to 
20kms _1 Vlsr range the mass is ~ 5 x 1O 4 M . Thus, the 
>1300Mq of extended gas traced by our broad- line NH3 
observations represents at least 5% of the total cloud mass. 



Cores 3 and 6: 

Cores 3 and 6 are not found towards any of the H.E.S.S. 
TeV gamma-ray sources. Their Vlsr values at ^-25kms -1 
are quite different to the other Cores in the region. The most 
likely connection is with the near-3-kpc spiral arm (with he- 
liocentric distance ^2 to 3 kpc), which has an expected value 
Vls r= —53.1 + 4.16Z [kms -1 ! for Galactic longitude / (see 
e.g. iDame et al. I [2008). Core 3 is possibly associated with 
the B8IV spectral type star HD 313632. Our new detec- 
tion of a H2O maser towards Core 3 may result from the 
envelope of HD 313632 or signal the presence of star forma- 
tion. From the Spitzer image in Figure a 24/mi feature is 
seen towards this core. For Core 6, the IR and radio sources 
IRAS 17555-2408 and PMN J1758-2405 are found - 2' and 
~ 4' distant. For both cores, our detection of NH3 could sig- 
nal some degree of star formation. Finally, Core 6 also lies 



just outside the TeV emission from HESS J1800-240C and 
thus is unlikely to be associated. For a 0.2 pc radius, both 
of these cores have masses ^600 to 900 M©, and densities 
-400 to 500 xlO 3 cm -3 . 



Cores 4 and 4a 

These cores are located towards the peak of the TeV 
gamma-ray source HESS J 1800- 240 A and are likely as- 
sociated with the Hll regions G6.225-0.569 (Core 4) an d 
G6.1-0.6 (Core 4a) (|Lockman II 19891 ; iKuchar fe Clarklll997h . 
Additional counterparts to Core 4a are the IR source, 
IRAS 17588-2358, also clearly visible in the Spit zer image 
(Figu re [7|, and a 1612 MHz OH maser ([Sevenster et al. I 
Il997h . Our detection of HC 3 N(3-2) towards Core 4a may 
also suggest it is at an earlier evolutionary stage than Core 4. 
Assuming a 0.2 pc radius we derived a mass and density of 
75 to 200 M and n H2 -45 to 100xl0 3 cm -3 for Core 4a 
and Core 4 respectively. Despite the reasonable amount 
of clumpy molecular gas traced by N an t en CO observa- 
tions (e.g. - 2.5 x 10 4 M from 12 CO (1-0) lAharonian et al. I 
2008b assuming a 2 kpc distance), we see no clear indica- 
tion for extended NH3 emission towards this region. We note 
however this region has only half the Mopra exposure com- 
pared to the other regions due to the lack of HOPS overlap. 

Triple Core and Core 5 
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The Triple Core represents the most complex of the re- 
gions we mapped, and comprises three NH3 peaks aligned 
in a SE to NW direction all of which are generally cen- 
tred on the TeV gamma-ray source HESS J1800-240B and 
the very IR-bright and energetic Hn region G5. 89-0. 39. Our 
12 mm observations are the largest scale mapping in dense 
molecular gas tracers so far of this enigmatic region. G5.89- 
0.39 actually comprises two a ctive star formation sites and 
following the nomenclature of iKim &; Kool ([2001) they are 
labeled G5.89-0.39A to the east, and G5.89-0.39B about 
2' to the west. The Triple Core SW NH3 core is associ- 
ated with Hllcore G5.89-0.39A, otherwise known as W28- 
A2 after its strong radio continuum emission. The ring- 
like features prominently visible in the Spitzer 8/rni emis- 
sion (see Figure [7J) are centred on G5.89-0.39A suggest- 
ing strong PAH molecular excitation from stellar pho- 
tons. The Triple Core Central NH3 core is associated with 
the UC-H11 region G5.89-0.3 9B from which st rong R76a 
RRL appears to be centred (|Kim &; Kooll200lh signalling 
strong ionisation of the surrounding molecular gas. The 
RRL H62a, H65a and H69a emission from our obser- 
vations also appears prominent towards G5.89-0.39B al- 
though in all three lines the emission is elongated towards 
G5.89-0.39A. The strongest H2O maser is also seen to- 
wards G5.89-0.39B. From the position-switched observa- 
tions, strong H2O maser emission with complex structures 
spanning very wide velocity coverage over 100 km s _1 are de- 
tected towards G5.89-0.39A and 60kms _1 in G5.89-0.39B. 
G5.89-0.39B is responsible for very energetic outflows and 



sec to arcmin scales, ( see e.g. lHarvev Forveille 



Churchwell et al.l Il990l: IGomezI Tl99ll: lAcord et al. 





1988 




1997 



Sollinsl 120041; [Thompson fe M acdonald Il999l; iKlaassenl 
2006; IKim fe Kooll200lL 120031 : iHunter et aTl 120081^) Prev i- 
ous sm a ll-scal e NH3 studies a r e disc u ssed by Gomez I (Il99lh ; 
I Wood I (|1993T ): lAcord et al.l (| 19971 ): IHunter et alT(|2008h . 
The Triple Core NW NH3 core is found a further 5' distant 
and may be linked to the M spectral type pulsating star 
V5561 Sgr or perhaps the natal gas from which this star was 
born. 

Core 5 appears to straddle the south east quadrant of 
the 8/mi IR shell or excitation ring of G5.89-0.39A and 
HESS J1800-240B, and is resolved into two components 
Core 5 NE and SW. Local peaks in CO emission overlapping 
Cor e 5 NE and SW are clearly visible (see F igure [6] and and 
also lLiszt et"aT1 (|2009h : IKim fe Kool (|2003h ). Core 5 NE is 
the coldest of the cores detected with Tk ~ 12 K. Here we de- 
tect only NH3 (1,1) emission in mapping, and only very weak 
(2,2) and (3,3) emission in the deep spectra. It is also one of 
the few sites where HC 5 N(10-9) is detected. In Core 5 SW 
we find Tk ~ 18 K and NH3 (2,2) and (3,3) emission being 
stronger than in the NE. There is also HCsN(3-2) emission 
when looking at the 5 to 20kms _1 velocity range in map- 
ping data. Overall, assuming a 0.2 pc core radius, we find our 
NH 3 (1,1) and (2,2) observations trace -100 to 450 M and 
density 71h 2 =80 to 300xl0 3 cm -3 for the individual cores 
in the Triple Core and Core 5 complexes. The core masses 
are in general agreement with their virial masses. 

Higher resoluti on NH 3 (3,3) observations of G5.39- 
0.39B by IGomezI (|l99lh with the VLA suggest a 
0.2 pc radius molecular envelope around G5.89-0.39B 
(Triple Core Central) tracing —30 M© assuming an abun- 



dance ratio Xnh 3 = 10 -6 . Using our abundance ratio this 
mass converts to — 1500 M©, about a factor 5 larger than 
our mass estimate using the NH3 (1,1) and (2,2) emission. 
This difference may arise from the fact that our analysis is 
not sensitive to the slightly broader (3,3) line. In addition 
(Purcell 2006) estimates a core mass (for radius 0.15 pc) us- 
ing SIMBA 1.2 mm continuum observations of 360 M© and 
derives similar virial masses from the linewidths of N2H + (1- 
0), H 13 CO (1-0), 13 CO (1-0) and CH 3 OH transitions, which 
are in general agreement with our results. 

Assuming the Triple Core and Core 5 complexes rep- 
resent extended sources, or at least the superposition of 
many unresolved point sources, we derive masses and den- 
sities (Table [5} of -3300 and 1300 M respectively with 
uh 2 ~ 0.7 x 10 3 cm -3 from the average NH 3 (1,1) and (2,2) 
emission. As for Core 2, such estimates do not consider 
the (3,3) emission seen towards Triple Core Central and SE 
and are therefore likely to underestimate the true extended 
mass. Nevertheless using these mass estimates the individual 
NH3 cores represent about 35% of the extended mass traced 
by NH3, for the two complexes, highlighting their gener- 
ally clumpy nature. We also find that the extended mass 
traced by NH3 represents only small fraction —5% of the 
total cloud mass, - 0.8 x 10 5 M , for the HESS J1800-240B 
region traced by the Nanten 12 C O (1-0) observations oy er 
the Vlsr=0 to 20kms _1 range ([Aharonian et al.ll2008bh . 

An interesting question is whether or not the W28 SNR 
shock has reached the southern molecular clouds, as it ap- 
pears to have done for Core 2 in the NE. Of interest therefore 
is the spatial distribution of any disruption in the molecu- 
lar clouds associated with the Triple Core and Core 5 com- 
plexes. The velocity dispersion image in Figure [9] clearly 
shows the broader NH3 gas is found concentrated towards 
the central star formation cores in constrast to Core 2 where 
the disruption appears to originate more from one side. As 
a result we would conclude there is no evidence (within the 
sensitivity limits of our mapping) for any external disruption 
of the southern molecular clouds. 



6 SUMMARY AND CONCLUSIONS 

We have used the Mopra 22m telescope for 12 mm line map- 
ping over a degree-scale area covering the dense molecu- 
lar gas towards the W28 SNR field. Our aim has been to 
probe the dense molecular cores of the gas spatially match- 
ing the four TeV gamma-ray source observed by the H.E.S.S. 
telescopes in this region. The wide 8 GHz bandpass of the 
Mopra telescopes allows to search for a wealth of molec- 
ular lines tracing dense gas. For our purpose, the emp- 
phasis has been on the inversion transitions of NH3 which 
allow relatively robust estimates of gas temperature and 
optical depth. Our observations combine data from dedi- 
cated scans and those from the the 12 mm HOPS project 
and reveal dense clumpy cores of NH3 towards most of the 
molecular cloud complexes overlapping the gamma-ray emis- 
sion. Additional 12 mm lines were detected, including the 
cyanopolyynes HC3N and HC5N, H2O masers, and radio 
recombination lines, which are prominent towards the UC- 
Hn complex G5.89-0.39. 

The NH3 cores are generally found in regions of 12 CO 
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peaks and mostly represent sites of star formation at various 
stages from pre-stellar cores to Hn regions. A standout ex- 
ception to this is the shocked molecular cloud on the north 
east boundary of the W28 SNR where we detect very broad 
Avi/2 > lOkms -1 NH3 emission up to the (9,9) transition 
with kinetic temperature >40K. The velocity dispersion of 
this broad gas further supports the idea that the SNR W28 
is the source of molecular disruption. 

Based on analysis of the NH3 (1,1) and (2,2) transitions 
we estimate total core masses for pure hydrogen in the range 
75 to 900 M© assuming a core radius of 0.2 pc. Molecular hy- 
drogen densities 71h 2 are in the range 40 to 500 xlO 3 cm -3 , 
however scaling factors for various core sizes included in Ta- 
ble [5] indicate that density varies with core radius while mass 
is more robust to changes in R. Obviously the estimates 
for core masses and densities are heavily dependent on the 
chosen NH3 abundance, however, the general agreement be- 
tween the derived mass and virial mass for the detected cores 
appears to vindicate our choice of Xnh 3 • Optical depths are 
in the range 1-4 and 0.2-1.2 for the (1,1) and (2,2) tran- 
sitions respectively. The north east NH3 core (Core 2) is 
extended and our radiative transfer analysis of the broad 
line components imply a total mass of at least 1300 M© and 
density 71h 2 > 0.7 x 10 3 cm -3 . This represents >5% of the 
total north east molecular mass traced by 12 CO (1-0) ob- 
servations. For the southern molecular clouds harbouring 
Hn regions, the velocity dispersion centres on the individual 
cores (Triple Core SE/Central/NW and Core 5 NE/SW) to- 
wards the centre of the molecular clouds. Thus at this stage 
we find no evidence (limited by the ^0.05 K per channel 
sensitivity of our mapping) for cloud disruption from exter- 
nal sources such as from W28 to the north, and thus it's 
likely the W28 shock has yet to reach the southern molecu- 
lar clouds. Of course this does not constrain the much faster 
transport of cosmic-rays from W28 to the southern clouds. 

Additionally we note the densities for Core 2, 
Triple Core and Core 5 (Table [5]) averaged over extended 
spherical regions are about a factor 10 below the critical 
density for NH3. This would suggest that much of the cores 
mass is probably contained within a smaller averaged vol- 
ume than that considered here, emphasising their clumpy 
nature. Another possible explanation for the lower densities 
could be the choice of spherical geometry, especially for the 
Core 2 region it is plausible that the shock has compressed 
the core deforming its geometry. 

This work is part of our ongoing study into the molecu- 
lar gas towards the W28 region, and help to further under- 
stand internal structures and dynamics of the molecular gas 
in the W28 SNR field. Deeper observations in 2010 reach- 
ing a ~0.02K channel sensitivity in the NH3 (1,1) to (6,6) 
(including (4,4)) transitions, have been made with Mopra 
in 12 mm towards the NE (Core 2). These deep observa- 
tions will permit pixel-by-pixel determination of gas param- 
eters such as density, temperature and line width. They will 
therefore permit a much more detailed probe of the effects 
of the SNR shock propagating into the Core 2 cloud, using 
theoretical predictions of the effects of shocks in molecu- 
lar clouds (IHollenbach fe McKeelfl979l : iDraine et all 1 19831 : 
iHollenbach &; McKeel ll989). Additionally, we recently com- 
pleted a survey in the 7 mm band to trace the disrupted 
and shocked gas with the SiO (1-0) and CS (1-0) lines. Fi- 
nally, we would also add that deeper gamma-ray observa- 



tions of the W28 sources will allow higher resolution gamma- 
ray imaging approaching the molecular core sizes revealed in 
this study. This, when combined with knowledge of molec- 
ular cloud structures on arcminute scales, will pave the way 
to probe the diffusion properties of cosmic-rays potentially 
producing the TeV gamma-ray emission. 
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